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ABSTRACT 

An understanding of the interaction between a moving conductor and a plasma is very impor- 
tant for the design of large structures in space, such as a space station or the tethered satellite 
system. It is well known that very large conducting objects which move slowly across mag- 
netic field lines radiate low frequency (Alfven) waves. In this paper we formulate the prob- 
lem in such a way that the radiation in all frequency bands can be computed. We then quote 
results of detailed calculations for spheres and cylinders of various sizes in the cold plasma 
approximation. In general, we find that in a plasma for which © 2 » © 2 and V 2 « c\ « c 2 , 
there is radiation in three frequency bands: CD < Q f , © M < © < Q e , co p < a < m uh , where V is 
the speed of the body, c A is the Alfven speed, c is the speed of light in vacuum, © is the fre- 
quency of the radiation, and Q„ © /A , Q e , co /; , and (O uh are the ion cyclotron, lower hybrid, elec- 
tron cyclotron, plasma, and upper hybrid frequencies, respectively. 

1. Statement of the problem 

We consider the problem of a conductor moving nonrelativistically through a magnetized 
plasma with velocity V = Ve r We call the rest frame of the plasma Z and the rest frame of 
the conductor Z . We assume that in Z far from the conductor the plasma is uniform and iso- 
tropic and has a magnetic field 


and an electric field 

Bq - B 0 e z 
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In Z' the fields far from the body 

are given by 


and 

B'o = So e z 

VxB 0 VB 0 

^ Gy 
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(3) 


(4) 


Inside the conductor, the electric field E' and the conduction current density J' c are related by 


t This paper is a shortened version of the paper of the same name that appeared in the Journal of Geophysical 
Research, Vol 91, No A9, pages 10, 117-10, 135, September 1, 1986. 
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In Z, Eq (5) takes the form 


J'c = CJE' 


(5) 


E = — 


VxB 0 

c 



( 6 ) 


We consider here the special case where V is perpendicular to B 0 . Our formulation can be 
easily modified to include the general case. The coordinate system that we use is shown in 
Fig. 1. 

A conductor placed in an electric field in vacuum becomes polarized so as to shield out 
the external field. The equipotential surfaces for a prolate spheroid embedded in an external 
electric field directed along its axis is shown in Fig. 2. The field has azimuthal symmetry 
about the axis of the spheroid. 

If the conductor is embedded in a plasma, the above solution cannot apply, since the 
component of the electric field parallel to the magnetic field will drive currents in the plasma. 
In particular, for an object is much longer in the direction of the external field than in the 
direction perpendicular to the field, such as the prolate spheroid shown in Fig. 2, the fields 
near the tips are much stronger than the field at infinity. 

Let L be a characteristic size of the body and c A the Alfven speed. If the body is large 
enough and is moving slowly enough (' V «c A and V/LQ, « 1), then the MHD approximation 
is applicable. The MHD model of the interaction between a moving conductor and a magnet- 
ized is based on the work of Drell, Foley and Ruderman (1965), who used it to explain the 
anomolously fast decay of the orbit of the Echo I weather satellite. The model did not gain 
wide acceptance until it was confirmed by in situ observations by Voyager I of the Alfven 
wing generated by the jovian satellite Io. Further theoretical work was done by Neubauer, 
and good agreement to the model was obtained in analysis of the magnetic field data by 
Acuna et al (1981) and of the plasma data by Barnett (1986). In this limit, the conductor 
radiates Alfven waves. Viewed in Z', there is a standing wave pattern consisting of two 
Alfven "wings" attached to the body. The wings have a cross-section determined by the 
shape of the body, and they extend in the direction of the Alfven characteristics, Vjjj, defined 
by 


v* = v± 


Bp 

V47ip 


(7) 


where p is the plasma mass density. A side view of the Alfven wing for a perfectly conduct- 
ing sphere is shown in Fig. 3. The electric field E' is zero inside the Alfven wings, while out- 
side them the electric field resembles the field that surrounds a conducting infinite cylinder 
emersed in a uniform electric field perpendicular to its axis. The external electric field is 
shielded out by charges on the surface of the wing. In addition, electric current flows along B 
on the surface of the wing. The current flows toward the sphere on the side of the wing that 
is negatively charged, crosses the magnetic field through the sphere, and flows away from the 
sphere on the side of the wing that is positively charged. A front view of the Alfven wing, 
showing schematically the charge and current density, is shown in Fig. 4. The electric field in 
a plane perpendicular to B through the center of the sphere is shown in Fig. 5. The plasma 
bulk velocity V' resembles the flow of an incompressible fluid around and infinite cylinder. 
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The streamlines in a plane perpendicular to B are shown in Fig. 6. An interesting property of 
the Alfvdn wing system is that a three dimensional problem has effectively become a two 
dimensional one. Instead of flowing around a spherical obstacle, the plasma must flow around 
a cylindrical one. 

Also shown in Figs. 3 and 4 is the direction of the Poynting vector S. Viewed in £, the 
conductor loses energy due to radiation. This radiation is analogous to classical Cherenkov 
radiation, which occurs whenever a particle moves through a dielectric faster than the phase 
speed of light in the dielectric. In the case of the Alfven wing, the conductor is moving 
accross the magnetic field faster than the Alfvdn waves can propagate in that direction (the 
phase velocity of Alfven waves perpendicular to B is zero). If the velocity of the conductor 
were greater than the Alfven speed, a shock would form instead of the Alfven wing. 

The question remains, are the MHD waves predicted by Drell et al the only waves gen- 
erated by a conductor moving through a plasma, or are there waves of higher frequencies? If 
other waves are possible, how does one compute the amount of energy radiated? 


2. Mathematical Formulation 


We now proceed to the formulation of the problem without the MHD frequency con- 
straint. We solve the problem using Fourier transforms. Our method incorporates the boun- 
dary conditions in an integral equation. We start by writing 


V 


x (V x E) + 


_ 1 _ j^E 

c 2 dP 


We now write J as a sum of two terms 


4k 9J 
c 2 dt 


( 8 ) 


J = J p (l-H) + J C H 
or 


(9) 


J «J C Jj 

where 


( 10 ) 


Js(x,0 = (J c (x,r) - J p (x,t))H(x,t) 
and 


( 11 ) 


1 1 inside the conductor 

0 elsewhere (12) 

J p is the current density in the plasma,while J c is the current density in the conductor. We 
now define the fourier transform /(k,co) of the function f(x,t) by 

/(k,co) = j j/(x,r)J = \ f(x,t ) e i( - k * ~ dt (13) 

J/>(k,co) is related to E(k,a>) through the dielectric tensor K by the relation 

J/>(k>0)) = —■ [K(k,co) - lj-E(k,co) (14) 

where I is the identity matrix, and K is the dielectric tensor of the plasma. Although the 
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computations that we quote at the end of this paper were performed using the cold plasma 
dielectric tensor, which is derived from the linearized equztions of motion, any plasma model 
can be used. Using (9) through (14), we can write the Fourier transform of (8) as 


T(k,co) • E(k,o)) = 


471/00 


Jc(k,co) 


(15) 


where T is the so-called dispersion tensor, defined by 

T = -k 2 1 + kk + — ^-K 

c 2 


(16) 


Eq. (16) can be inverted to solve for E in terms of J s 

E(k,co) = - -^T^M-J^k,®) 
cr 


(17) 


Once E(k,co) is known, the magnetic field perturbation B, the plasma current density J p , and 
the electric charge density p can easily be found from the Fourier transforms of Maxwell’s 
equations. We seek an equation for J s . To derive it, we multiply Eq. (6) by H(x,t ) and use 
Eq. (11) to eliminate J c . The result is in the linlearized approximation 
V x B = V x B 0 is 




VxBq' 


Js(x,r) = 

a 

E(x,r) + 

c 

k. J 

- J pM + p c (x,r)V 




(18) 


where we have linearized by writing B 0 instead of B. We note that the Fourier transforms of 
the fields on the right side of Eq. (18) can be written in terms of the Fourier transform of J r 
The result is the desired integral equation for J s . 



(k,co) 


r = 0 


for x in the plasma 


(19a) 



VxB 0 

for x in the conductor 

c 


(19b) 


where the tensor W is defined by 


W = 


_ 4tu'co t _ 1 _1_ 


(k 2 - -^-)I - kk 


r-1 


( 20 ) 


We are led to an integral equation by the facts that the relation between E and J 5 is an alge- 
braic relation in (k,co) space, while the location of the conductor, as expressed by H(x,t), is 
simple in ordinary (x,r) space. Notice that the boundary conditions are automatically included 
in our formulation. 


3. Properties of the solutions 

The solution to (19) and (20) contains a complete description of the fields surrounding a 
conductor moving through a magnetized plasma. Unfortunately, this integral equation is 
difficult to solve. To obtain numerical estimates, we therefore take the alternate approach of 
trying to study the properties of the solutions. If we assume that we know the source current 
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Js> the electric field is given by the Fourier transform of (17). When one uses a complex con- 
tour integration to evaluate the integral, the only contributions come from the poles of the 
integrand, which occur at the zeros of the equation 

det T 1 = 0 ( 2 i) 

Equation (21) is the familiar dispersion relation for the plasma. In the present case, we seek a 
solution which is independent of time in 2', which implies that in 2, 

CO = k-V (22) 

For the geometry described in (1) - (4), (22) takes the form 

co = k x V ( 23 ) 

Furthermore, since we are primarily interested in the power radiated, we need to study the 
behavior of the solutions far from the origin. In particular, we are interested only in those 
modes for which k and co are both real, since only these modes can transport energy to 
infinity. In low Earth equatorial orbit, the plasma and the orbital velocity obey the following 
inequalities ° 

<4 » (24) 

and 

V 2 « c\ « c 2 (25) 

For such conditions, (21) and (23) have solutions for real k and co in only three frequency 
bands 

03 < Q Band I 

co_;jj < co < Q. e Band II (26) 

CO), < co < co„ A Band III 

Band I is the MHD band discussed by Drell et al; bands II and III are new. This radiation is 
analogous to classical Cherenkov radiation. We have already pointed out that a conductor 
moving through a magnetized plasma only radiates plasma waves whose phase speed in some 
direction is slower than the speed of the body. For a body that moves slowly compared to the 
Alfven speed, this occurs only in the frequency bands described by (26). A polar plot of the 
phase speed versus the angle between k and B 0 for the mode that is excited in each of these 
three bands is shown in Figs. 7-9. Note that the phase velocity vanishes in some direction for 
each of these three cases. 


4. Calculation of radiated power 

Having identified the frequency bands within which power is radiated, we now have to 
estimate the amount of power that will be radiated into each band. To do this, we guess the 
spatial dependence of the source current distribution and use Poynting’s theorem to compute 
its magnitude. Poynting’s theorem can be written as follows 


Prod = -W- ~ 


*Vem 

dt 


(27) 
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where 


rad 


c 

4tc 

J E x B • da 

(28) 

w = j 

E Jdx 

(29) 

1 

8rt 

j(E 2 + B 2 ) dx 

(30) 


where P rad is the power radiated, U E u is the electro-magnetic energy, and W is the mechani- 
cal work done by the fields. We choose for our volume of integration a large rectangular 
prism with the conductor at its center, and we perform all of the calculations in frame Z. If 

bU EM 

we make the box large enough, the term — - — can be neglected, due to the fact that the 

at 

solution we seek is independent of time in frame Z'. To evaluate the remaining integrals, we 
use the following ansatz; assume that the source current can be written as 


J s (x) = - 1 f(x) e y (31) 

where /(x) is a known function, the form of which one guesses from the geometry of the con- 
ductor. Now consider (28). Since both E and B are linear functions of J s , P rad is proportional 
to I 2 . To evaluate (29), we use Ohm’s law (6) for E. Since J c is proportional to J 5 , (29) 
results in two terms, one proportional to I and one proportional to I 2 . We can therefore write 
(27)-(29) as 


I 2 Z rad = /(<)>- I 2 R) (32) 

where <j> has dimensions of electric potential, and we have assumed that J s = J c . Z rad , R, and <}> 
can be computed using (31), (17), (13), (28) and (29). Equation (32) can be interpreted as an 
electrical circuit analog, where <{> is the voltage, I is the current, and Z rad and R are the 
resistences of two resistors in series. The problem of computing the radiated power thus 
reduces to evaluating Z rad , given by 


with E and B given by 


and 


z «* = ~hi \ < E x B >- * 


E(x,r) = J " 1 


47tt(0 T _1 


/(k)e. 



(33) 


(34) 


(35) 


where /( k) is the Fourier transform of/(x) defined in (31). We have evaluated (33) explicitly 
for two different geometries. The first case is a sphere of radius a. The second is a cylinder 
of radius b and length L whose axis points in the e^, direction. In both cases we have assumed 
a current density of the form given in (31), with/(x) a constant inside the object and zero out- 
side it. In both cases Z rad can be expressed as a sum of three terms 
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^rad ~ 2/ + Z 77 + Zjjj (36) 

where Z 7 , Z 77 , and Z n are the contributions from bands I, II, and III, respectively. Analytic 
expressions for Z as a function of a and of b and L can be found for several limiting cases. 
These cases are summarized in Table 1. The results of numerical calculations of Z 7 , Z 77 , Z Ib 
an< * P ra d are shown in Figs. 10-13, and the plasma parameters used in the calculations are 
given in Table 2. Of considerable interest are the results for a long, thin wire, since this might 
have applications to tethered satellite systems. If one wants to estimate the power that can be 
generated using a TSS, one can simply consider R to be the load resistence. Due to the con- 
tributions from band H, the impedence is much higher than expected in this case. A cylinder 
with L = 10 km and b - lcm has a radiation impedence of nearly 10 5 ohms! These results 
suggest that previous power estimates for passive TSS systems might be much too high. One 
must use care when applying these results to the TSS, however. In particular, effects due to 
the size of the subsatellite and local plasma clouds from "plasma contactors" can be mocked 
up by using an "effective" b. 

Our formulation can also be used to estimate the radiated power for active systems. One 
must then consider / in (31) as a known source strength Some calculations of that sort are dis- 
cussed in the article by Hastings et al. 

5. Summary 

We have considered the problem of the interaction between a moving conductor and a 
magnetized plasma. We have shown that steady state solutions exist for which the body is 
surrounded by a system of standing waves. Wave modes are excited for which the phase 
speed in some direction is less ’than the speed of the conductor. The process is analogous to 
Cherenkov radiation, or to the formation of a shock wave surrounding a body in supersonic 
motion through a gas. We have extimated the power radiated for spherical and cylindrical 
bodies moving perpendicular to the background magnetic field. Our results suggest that, due 
to the high radiation impedence in band n, a passive electrodynamic tethered satellite system 
will not draw large currents. 
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TABLE 1. Radiation Impedence of Conductors Moving Through a 
Plasma Perpendicular to the Magnetic Field 


Size 


Impedence 


Sphere of Radius a 


Region I 


V 

a » — 
Hi 

V 

V 

a » ■ 


Region II 


Vlh 

V V 

— « a « 

V 


Z -~ = Z 
c z = z o 

7 -^li7 

Z,_ l5F Z ° 

Z - *« Vc * 7 
Z " 3 Z ° 

Z - ^ :tr< z 

Z - — Z 

Z " VF Z ° 


Reaion III 


a » 


iO LH C 


Zut — 


AVc, 


15a*a> UJ ' 1 Z ° 


Cylinder of Length L and Radius b With Axis Perpendicular 


to Both V and B„ 


Reaion I 


Reaion II 


Region III 


. V V 
b « — L « — 

A, 

b « — L» — 

G| 

L» — b « — 


_ jrLflj 

z ‘ = 7F z ° 


z, 

z, 


.[,„(^ + 0 , 7 ] Zo 

»=M ,n fe) + U16 ) 


L » — £> » — Z„, = 




tu„ 


cu. 


cibcouj) 


J 0 


We have expressed the impedence in terms of Z 0 = 2c Jc 2 in Gaus- 
sian units. For conditions typical of low earth orbit, Z 0 = 6.7 
x 10"‘* s/cm = 0.06£1 If rj = = 1/172 (oxygen plasma). 

c< = 300 km/s and B 0 = 0.33 G. 


TABLE 2. Plasma Parameters Used in.Numerical Calculations 


Parameter Value 


B 

0.33 G 

”, 

3.6 x I0 J cm" J 

1 1 

1/172 

K 

7.3 km/s 

* A 

300 km/s 

M A 

0.024 


3.4 x 10 7 s*‘ 

0. 

5.9 x 10 6 s*‘ 

Um 

3.5 x IO^s" 1 


2.0 x I0 J s"‘ 

UK 

1.2 mm 

h'^LH 

21 cm 

UK 

36 m 
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Fig. 4. Alfven wing for perfectly 

conducting sphere-front view 



5. Electric field lines in Fig. 6. Stream lines in equatorial plane for 

equatorial plane for flow around perfectly conducting sphere 

Derfectly conducting sphere 
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Fig. 10. Impedence vs. radius for 
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